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ABSTRACT 

We  performed  simulations  of  initial  stages  of  a 
carbon  nanotube  growth  catalyzed  by  iron  particles  using 
all-electron  density  functional  theory  with  generalized 
gradient  approximation  for  the  exchange- correlation 
functional.  The  systems  studied  are  Fe4  -i-  C^  for  n=7-25, 
Fcio  +  Ci5,  and  Fe^oCis  -i-  IOC.  In  addition,  we  performed 
a  detailed  study  of  development  of  different  isomers  in 
the  Cn  series,  n=7-25. 

1.  INTRODUCTION 

Single- walled  carbon  nanotubes  (SWCNT)  are 
expected  to  have  an  enormous  impact  on  various 
technological  areas  related  to  fabrication  of  composite 
materials,  sensors,  hydrogen  storages,  protective  coatings, 
and  computer  memories.  The  SWCNTs  are  produced 
typically  from  hydrocarbons  [1],  alcohol  [2],  and  graphite 
[3].  The  method  of  catalytic  chemical  vapor  deposition 
(CCVD)  [1]  makes  use  of  high  temperatures  to  produce 
carbon  and  a  supported  catalyst  to  initiate  the  carbon 
nanotube  growth.  The  catalyst  needs  restoration  after  each 
cycle.  Carbon  nanotubes  grown  using  this  method  are 
often  multi- walled,  while  single- walled  nanotubes  are  of 
more  interest  from  a  technological  point  of  view.  On  the 
contrary,  the  high-pressure  high-temperature  HiPco  [4] 
method  is  continuous  and  its  floating  or  gas-phase  catalyst 
is  formed  in  situ  from  iron  pentacarbonyl  Fe(CO)5.  Using 
carbon  monoxide  as  feedstock  is  rather  intriguing,  because 
CO  is  one  of  the  most  stable  diatomic s  with  the  measured 
bond  strength  of  11.09  eV.  Carbon  nanotubes  produced  in 
this  process  are  always  single-walled  and  are  believed  to 
nucleate  and  grow  via  the  Boudouard  disproportionation 
reaction  CO  -H  CO  ^  Cswcnt  +  CO2.  However,  the 
mechanism  of  carbon  nanotube  growth  is  not  well 
understood. 

In  order  to  gain  insight  into  carbon  nucleation  in 
the  HiPco  process,  we  have  performed  all-electron  density 
functional  theory  calculations  with  a  generalized  gradient 
approximation  (DFT-GGA)  on  [5],  FCnCO  (n<6)  [6], 
and  Fe4Cn(CO)jn  (n-i-m<3,  both  C  and  CO  are  attached  to 
FQ4)  [7]  clusters,  and  estimated  the  energies  of  Boudouard 
CO  disproportionation  reactions  Fe4Cn(CO)ni+ 
CO^Fe4Cn+i(CO)ni-i+  CO2.  The  energies  found  are 
relatively  small  and  are  in  the  range  from  -0.26  eV  to 
-1-0.74  eV.  Next,  we  explored  [8]  the  thermodynamics  of 
CO  attachment  to  a  carbon  atom  that  had  already 


precipitated  on  an  iron  cluster  due  to  a  preceding 
Boudouard  reaction  and  evaluated  the  energetics  of  CO 
disproportionation  reactions  Fe4Cn(CO)in  +  CO  ^ 
Fe4Cn+i(CO)in.i  +  CO2  for  higher  coverage  (n-i-m<5).  We 
choose  Fe4',  Fe4,  and  Fe4'^  clusters  since  the  Fe4'^  cluster  is 
known  to  catalyze  the  growth  of  benzene  from  ethylene 
and  cyclopropane  in  a  low-pressure  gas-phase  process. 
The  number  of  iron  atoms  in  the  clusters  formed  in  the 
Hipco  process  ranges  from  '-lO  to  '-300  and  it  is  not  clear 
if  smaller  iron  clusters  Fe3  to  Fciq  can  technologically  be 
effective  for  the  SWCNT  growth  because  of  different 
restrictions  such  as  coalescence  of  small  clusters. 
However,  one  can  anticipate  that  our  computations  using 
computationally  less  demanding  Fe4  clusters  are  capable 
of  reproducing  the  essentials  in-initiating  the  SWCNT 
growth.  According  to  the  results  of  our  previous 
computations,  dimerization  of  carbon  atoms  is  preferable 
in  Fe4C4,  Fe4C5  and  its  ions,  while  C3  trimers  form  in  the 
ground  states  of  Fe4C6  and  its  ions.  Since  attachment  of 
COs  followed  by  subsequent  Boudouard  reaction  leads  to 
precipitation  of  carbon  atoms,  we  considered  next  the 
trends  in  rearrangement  of  carbon  atoms  in  the  lowest 
energy  isomers  of  clusters  Fe4Cn,  n=7-20  [9].  This  work 
discusses  some  results  from  Ref.  [9]  and  presents  the 
results  of  our  simulations  using  and  Fe4Ci5  clusters. 

2.  COMPUTATIONAL  DETAILS 

The  Gaussian  98  and  03  programs  [10,11]  were  used. 
We  have  used  the  6-311-t-G*  basis  set 
(15sllp6dlf)/[10s7p4dlf]  [12]  and  (12s6pld)/[5s4pld] 
[13]  for  Fe  and  C,  respectively.  In  some  computations 
involving  Fcio  and  Fe4Ci5,  we  used  a  smaller  6-3 IG  basis 
for  carbon  atoms.  Our  previous  study  [14-16]  of  bare  iron 
clusters  and  3d-metal  dimers  showed  that  results  obtained 
using  many  of  the  DFT-GGA  methods  included  in 
Gaussian  98  and  03  programs  are  rather  similar;  however, 
the  BPW91  vibrational  frequencies  appear  to  be  less 
sensitive  to  the  quality  of  the  grid  used  in  the  numerical 
integration  than  some  of  the  other  functionals.  On  this 
ground,  we  choose  the  BPW91  method,  where  the 
exchange- correlation  functional  is  comprised  of  the 
Becke's  exchange  [17]  and  Perdew- Wang's  correlation 
[18]. 

The  geometry  of  each  cluster  was  optimized 
without  imposing  any  symmetry  constraints.  Iron  clusters 
possess  high  spin  multiplicities  and  our  optimizations  of 
the  clusters  interacting  with  carbon  species  were 
performed  in  the  range  of  spin  multiplicities  enclosing  the 
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ground- state  spin  multiplicity  of  the  corresponding  bare 
iron  cluster.  Each  geometry  optimization  was  followed  by 
the  calculation  of  the  harmonic  vibrational  frequencies 
using  analytical  second  derivatives,  in  order  to  confirm 
that  the  optimized  geometry  corresponds  to  a  minimum. 
We  computed  atomic  spin  densities  using  Mulliken  [19] 
approach. 


3.  GEOMETRICAL  STRUCTURES  OF 
Fe4C„ 


Our  optimizations  revealed  a  large  number  of 
isomers  for  each  Fe4Cn  with  a  general  trend  of  the  carbon 
ring  formation  The  rings  may  contain  also  iron  atoms  in 
Fe4C6-Fe4Ci5.  The  electronic  states  of  C^  possessing  ring- 
type  geometries  in  the  gas  phase  have  the  lowest  total 
energies  [20]  beginning  with  Ciq.  We  found  that  the 
formation  of  pure  carbon  rings  in  the  lowest  total  energy 
of  the  Fe4Cn  clusters  begins  with  n=15-16.  Figure  1 
presents  the  results  of  our  computations  for  clusters  Fe4C8 
and  Fe4C9  and  shows  free  carbon  clusters  Cg  and  C9 
optimized  at  the  same  level  of  theory. 


FG4C9  Cg 


2S+1=15  ground  state  ^Sigma^  0.0  eV 


Fig.  1.  Top  and  side  views  of  the  ground- state  Fe4  clusters 
with  attached  carbon  Cg  and  C9  species.  The  ground  and 
first  excited  ring  states  of  the  latter  are  shown  in  the  right- 
hand  side.  Bond  lengths  are  in  A,  excess  electronic 
densities  at  iron  sites  are  in  electrons. 


As  is  seen,  carbon  atoms  form  ring  structures 
containing  two  iron  atoms  in  Fe4C8  and  in  Fe4C9,  while 
the  ground  states  of  Cg  and  C9  are  linear.  The  ring 
structures  of  Cg  and  C9  are  above  the  ground  states  by 
about  0.6  eV.  The  ground- state  of  the  Cg'  and  C9'  anions 
are  also  linear.  As  for  the  C^'  anions,  there  are  near- linear 


states  that  were  a  subject  of  experimental  [21]  and 
theoretical  [22]  studies.  According  to  Ref.  [22]  and  our 
computations  performed  with  a  larger  basis  set  than  that 
used  in  Ref.  [22],  the  anion  ring  structure  is  more  stable 
than  a  nearly  linear  one  up  to  n  =  14.  (Note  that 
optimizations  of  C^'  anions  imposing  Dooh  constraints 
produce  transition  states  beginning  with  n=10.)  It  appears 
that  this  number  is  critical  for  formation  of  carbon  ring 
structures  on  the  top  of  an  iron  catalyst  since  there  is  a 
charge  transfer  from  iron  particles  to  a  more 
electronegative  carbon  chain.  Figure  2  shows  that  Fe4Cio 
contains  a  ring  with  two  inserted  iron  atoms  while  Fe4Ci4 
forms  a  ring  with  one  inserted  iron 
atom. 


^10 


2S+1  =  11  ground  state  C2h  Near-linear  2Ay-2.02eV 


^14 


2S+1  =  11  ground  state  ground  state 

Fig.  2.  Top  and  side  views  of  the  ground-state  Fe4Cio  and 
Fe4Ci4  clusters  together  with  the  ground  states  of  of  C^q, 


C]^Q ,  and  C]^4. 


In  Fe4Ci5,  there  is  a  competition  between 
configurations  containing  a  pure  C15  ring  and  a  C15  ring 
with  an  iron  atom  inserted  inside  the  ring.  The  pure 
carbon  ring  structure  on  the  top  of  Fe4  becomes  favorable 
in  Fe4  -H  The  lowest  energy  state  of  Fe4  -1-  C20 
correspond  to  a  C19  ring  around  the  cluster  and  one  carbon 
atom  attached  to  a  face  of  the  iron  cluster. 

In  order  to  gain  insight  which  geometrical 
structure  corresponds  to  the  lowest  energy  state  in  Q  we 
performed  optimizations  [9]  of  Q  up  to  n=25.  The  first 
fullerene-type  structure  with  12  carbon  pentagons  appears 
at  n=20,  but  the  ring  structure  is  the  lowest  energy  one, 
followed  by  a  bowl  (flake-type)  structure,  and  a  fullerene 
(or  cage)  structure.  Increasing  the  n  number  to  n=25,  we 
found  the  ring  structures  to  correspond  to  the  lowest 
energy  states  for  n  =  21,  22,  23,  and  25,  while  the  lowest 
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energy  state  corresponds  to  a  flake  composed  of  seven 
hexagons  for  n=24.  At  this  size,  there  appears  the  smallest 
fullerene  possessing  two  hexagons  in  addition  to  12 
pentagons.  The  results  of  our  optimizations  for  n=24  are 
presented  in  Fig.  4.  There  are  a  large  number  of  different 
shape  isomers;  some  of  them  are  shown  in  the  figure.  The 
triplet  states  are  higher  for  each  flake,  ring,  or  fullerene 
structures. 


C 


16 


2s+l  =  13  ground  state 


^Ag  ground  state 


2s+l  =11  ground  state 


^Ag  ground  State 


Fig.  3.  Top  and  side  views  of  the  ground- state  Fe4  clusters 
interacting  with  and  C20. 

4.  ISOMERS  OF  Fcjo  AND  FcjoCis 

To  explore  the  influence  of  an  iron  cluster  size, 
we  have  optimized  several  isomers  of  Fciq.  Next,  we 
added  a  C15  ring  to  one  of  the  isomers  in  different 
positions:  on  the  top  of  long  and  short  sides  and  wrapped 
around  the  center  part  of  the  cluster.  The  results  of  our 
optimizations  are  presented  in  Fig.  5. 

We  found  that  the  carbon  ring  placed  along  the 
larger  side  has  acquired  one  inserted  iron  atom.  This 
configuration  corresponds  to  the  lowest  energy  state.  The 
ring  placed  over  a  shorter  side  of  the  iron  cluster  slightly 
changes  the  perfect  ring  shape  and  attaches  perpendicular 
to  the  iron  cluster.  The  C15  ring  placed  around  the  center 
shifts  during  lengthy  optimization  cycles  to  the  cluster  tip 
while  the  initial  bi-octahedral  Fe^o  configuration 


transforms  to  another  isomer  configuration.  The  carbon 
ring  placed  above  the  tip  is  tilted  to  be  nearly  parallel  to  a 
Fe3  edge.  The  latter  two  Fe^oCis  isomers  are  nearly 
degenerate  in  total  energy  and  are  above  the  ground  state 
by  about  3  eV.  The  binding  energy  gain  in  the  ground 
state  apparently  comes  from  breaking  a  C-C  bond  and  the 
formation  of  two  Fe-C  bonds.  This  is  apparently  related 
with  the  shape  of  the  surface  and  does  not  happen  when 
Ci5  is  attached  to  a  shorter  side 

Isomers  of  C24 


Distorted  ring 
U(Ci)  +3.10  eV 


Flake 

lAg(C2h)  0.0  eV 


Cap 

1a  (C2)  +6.88  eV 


Non-interacting  C,2rings 
1a(Ci)  +7.05  eV 


Ring 

^Ag  (C2h)  +0.28  eV 


Fullerene 
1a  (C2)  +0,45  eV 


Fig.  4.  .Different  isomers  of  C24. 


One  could  speculate  that  the  upward  growth  on 
the  tilted  ring  might  produce  a  carbon  nanotube  with  a 
certain  chirality.  That  is,  the  part  of  the  cluster  surface, 
which  initiates  the  carbon  ring,  is  related  to  the  chirality 
of  a  SWNT  grown  from  the  corresponding  carbon  ring. 

5.  SIMULATED  +  IOC  GROWTH 


In  order  to  gain  an  insight  how  carbon  atoms 
attach  to  a  ring  formed  on  an  iron  particle,  we  chose  a 
F^4Ci5  isomer  with  a  ring  C15  pattern  and  added  10  carbon 
atoms  randomly  placed  at  the  distance  of  5  A  above  the 
ring.  Initial  and  final  configurations  together  with  several 
intermediate  configurations  obtained  at  different 
optimization  steps  are  presented  in  Fig.  6. 

It  is  found  that  carbon  atoms  may  form  first  dimers  and 
trimers  in  correspondence  with  experimental  observations 
of  such  a  formation  in  arc  discharges.  Next,  chains  are 
formed,  while  the  final  configuration  correspond  to  a  cage 
with  an  open  15-carbon  ring. 

CONCLUSION 


Our  simulations  at  a  rather  reliable  DFT-GGA  level  show 
that  the  prime  role  of  a  catalyst  is  to  assemble  carbon 
atoms  in  the  optimal  gas-phase  configuration  of  C^  that 
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corresponds  to  a  ring  configuration  for  a  rather  large  n. 
The  surface  topology  of  a  catalyst  particle  is  most  likely 


TOP  VIEW 


SIDE  VIEW 


M=29  0.0  eV 


M=29  +2.81  eV 


M=29  0.0  eV 


TOP  VIEW 


SIDE  VIEW 


M=31  +0.01  eV 


M=31  +2.89  eV 


M=31  +0.46  eV 


Fig.  5.  Different  lowest  energy  isomers  of  Fe^o  Fe^oCi 


M=ll  0.0  eV 

Q 


M=ll  +7.6  eV 


M=ll  +13.8eV 


M=1 1  +44.2  eV 


M=ll  +34.0  eV 


M=ll  +23.1  eV 


responsible  for  the  chirality  of  a  carbon  nanotube  grown 
from  a  seed  ring. 
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